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1. BIG-BANG NUCLEOSYNTHESIS 

Revised October 2005 by B.D. Fields (Univ. of Illinois) and S. Sarkar (Univ. of Oxford). 

Big-bang nucleosynthesis (BBN) offers the deepest reliable probe of the early universe, 
being based on well-understood Standard Model physics [1-4]. Predictions of the 
abundances of the light elements, D, 3 He, 4 He, and 7 Li, synthesized at the end of the 
"first three minutes" are in good overall agreement with the primordial abundances 
inferred from observational data, thus validating the standard hot big-bang cosmology 
(see [5] for a review). This is particularly impressive given that these abundances 
span nine orders of magnitude — from 4 He/H ~ 0.08 down to 7 Li/H ~ 10 -10 (ratios 
by number). Thus BBN provides powerful constraints on possible deviations from the 
standard cosmology [2], and on new physics beyond the Standard Model [3]. 

1.1. Theory 

The synthesis of the light elements is sensitive to physical conditions in the early 
radiation-dominated era at temperatures T ^ 1 MeV, corresponding to an age t ^ 1 s. 
At higher temperatures, weak interactions were in thermal equilibrium, thus fixing 
the ratio of the neutron and proton number densities to be n/p = e~Q/ T , where 
Q = 1.293 MeV is the neutron- proton mass difference. As the temperature dropped, 
the neutron- proton inter-conversion rate, V^^p ~ GpT 5 , fell faster than the Hubble 
expansion rate, H ~ V 9*&N T 2 , where g* counts the number of relativistic particle 
species determining the energy density in radiation. This resulted in departure from 
chemical equilibrium ("freeze-out") at Tf r ~ (fi^GWGj^) 1 / 6 ~ 1 MeV. The neutron 
fraction at this time, n/p = e~Q/ T h ~ 1/6, is thus sensitive to every known physical 
interaction, since Q is determined by both strong and electromagnetic interactions while 
Tf r depends on the weak as well as gravitational interactions. Moreover the sensitivity 
to the Hubble expansion rate affords a probe of e.g. the number of relativistic neutrino 
species [6]. After freeze-out the neutrons were free to /3-decay so the neutron fraction 
dropped to ~ 1/7 by the time nuclear reactions began. A simplified analytic model of 
freeze-out yields the n/p ratio to an acuracy of ~ 1% [7,8]. 

The rates of these reactions depend on the density of baryons (strictly speaking, 
nucleons), which is usually expressed normalized to the relic blackbody photon density 
as n = ng/n-y. As we shall see, all the light-element abundances can be explained with 
7710 = i] x 10 10 in the range 4.7-6.5 (95% CL). With n 7 fixed by the present CMB 
temperature 2.725 K (see Cosmic Microwave Background review), this can be stated as the 
allowed range for the baryon mass density today, pq = (3.2-4.5) x 10 -31 gem -3 , or as the 
baryonic fraction of the critical density, fig = ps/Pcrit — r ?l0^ _2 /274 = (0.017-0.024) /i -2 , 
where h = i^o/100 kms Mpc -1 = 0.72 ± 0.08 is the present Hubble parameter (see 
Cosmological Parameters review). 

The nucleosynthesis chain begins with the formation of deuterium in the process 
p(n, 7)D. However, photo-dissociation by the high number density of photons delays 
production of deuterium (and other complex nuclei) well after T drops below the binding 
energy of deuterium, = 2.23 MeV. The quantity r/ _1 e _ziD / r , i.e. the number of 
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photons per baryon above the deuterium photo-dissociation threshold, falls below unity at 
T ~ 0.1 MeV; nuclei can then begin to form without being immediately photo-dissociated 
again. Only 2-body reactions such as D(p, 7) 3 He, 3 He(D,p) 4 He, are important because 
the density has become rather low by this time. 

Nearly all the surviving neutrons when nucleosynthesis begins end up bound in the 
most stable light element 4 He. Heavier nuclei do not form in any significant quantity 
both because of the absence of stable nuclei with mass number 5 or 8 (which impedes 
nucleosynthesis via n 4 He, p 4 He or 4 He 4 He reactions) and the large Coulomb barriers for 
reactions such as T( 4 He, 7) 7 Li and 3 He( 4 He, 7) 7 Be. Hence the primordial mass fraction of 
4 He, conventionally referred to as Y p , can be estimated by the simple counting argument 

Y p = ^Uo.25. (1.1) 
p 1 + n/p V ' 

There is little sensitivity here to the actual nuclear reaction rates, which are however 
important in determining the other "left-over" abundances: D and He at the level of a 
few times 10 -5 by number relative to H, and 7 Li/H at the level of about 10 -10 (when 7710 
is in the range 1-10). These values can be understood in terms of approximate analytic 
arguments [8,9]. The experimental parameter most important in determining Y p is the 
neutron lifetime, r n , which normalizes (the inverse of) Tr^p. (This is not fully determined 
by Gp alone since neutrons and protons also have strong interactions, the effects of which 
cannot be calculated very precisely.) The experimental uncertainty in r n used to be a 
source of concern but has recently been reduced substantially: r n = 885.7 ± 0.8 s (see N 
Baryons listing). 

The elemental abundances, calculated using the (publicly available [10]) Wagoner 
code [1,11], are shown in Fig. 1.1 as a function of 7710- The 4 He curve includes small 
corrections due to radiative processes at zero and finite temperature [12], non-equilibrium 
neutrino heating during annihilation [13], and finite nucleon mass effects [14]; the 
range reflects primarily the lcr uncertainty in the neutron lifetime. The spread in the 
curves for D, 3 He and 7 Li corresponds to the la uncertainties in nuclear cross sections 
estimated by Monte Carlo methods [15-16]. Recently the input nuclear data have been 
carefully reassessed [17-21], leading to improved precision in the abundance predictions. 
Polynomial fits to the predicted abundances and the error correlation matrix have been 
given [16,22]. The boxes in Fig. 1.1 show the observationally inferred primordial 
abundances with their associated statistical and systematic uncertainties, as discussed 
below. 



1.2. Light Element Abundances 

BBN theory predicts the universal abundances of D, 3 He, 4 He, and 7 Li, which are 
essentially determined by t ~ 180 s. Abundances are however observed at much later 
epochs, after stellar nucleosynthesis has commenced. The ejected remains of this stellar 
processing can alter the light element abundances from their primordial values, but also 
produce heavy elements such as C, N, O, and Fe ("metals"). Thus one seeks astrophysical 
sites with low metal abundances, in order to measure light element abundances which 
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Figure 1.1: The abundances of He, D, He and Li as predicted by the standard 
model of big-bang nucleosynthesis. Boxes indicate the observed light element 
abundances (smaller boxes: 2a statistical errors; larger boxes: ±2cx statistical and 
systematic errors). The narrow vertical band indicates the CMB measure of the 
cosmic baryon density. See full-color version on color pages at end of book. 

are closer to primordial. For all of the light elements, systematic errors are an important 
and often dominant limitation to the precision with which primordial abundances can be 
inferred. 
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In recent years, high-resolution spectra have revealed the presence of D in high-redshift, 
low-metallicity quasar absorption systems (QAS), via its isotope-shifted Lyman-a 
absorption [23-27]. It is believed that there are no astrophysical sources of deuterium 
[28], so any detection provides a lower limit to primordial D/H and thus an upper limit 
on 77; for example, the local interstellar value of D/H = (1.5 ± 0.1) x 10 -5 [29] requires 
7710 < 9. Recent observations find an unexpected scatter of a factor of ~ 2 [30] so 
interstellar D may have been depleted by stellar processing. However, for the high-redshift 
systems, conventional models of galactic nucleosynthesis (chemical evolution) do not 
predict D/H depletion [31]. 

The observed extragalactic D values are bracketed by the non-detection of D in a 
high-redshift system, D/H < 6.7 x 10 -5 at la [32], and low values in some (damped 
Lyman-a) systems [24,25]. Averaging the 5 most precise observations of deuterium 
in QAS gives D/H = (2.78 ± 0.15) x 10 -5 , where the error is statistical only [23,26]. 
However there remains concern over systematic errors, the dispersion between the values 
being much larger than is expected from the individual measurement errors (x 2 = 16.5 
for 7y = 4 d.o.f.). Increasing the error by a factor a/x 2 /^ gives, as shown on Fig. 1.1: 



We observe He in clouds of ionized hydrogen (H II regions), the most metal-poor 
of which are in dwarf galaxies. There is now a large body of data on 4 He and CNO in 
these systems [33]. These data confirm that the small stellar contribution to helium is 
positively correlated with metal production. Extrapolating to zero metallicity gives the 
primordial 4 He abundance [34] 



Here the latter error is a careful (and significantly enlarged) estimate of the systematic 
uncertainties; these dominate, and is based on the scatter in different analyses of the 
physical properties of the H II regions [33,34]. Other extrapolations to zero metallicity 
give Y p = 0.2443 ± 0.0015 [33], and Y p = 0.2391 ± 0.0020 [35]. These are consistent 
(given the systematic errors) with the above estimate [34], which appears in Fig. 1.1. 

The systems best suited for Li observations are metal-poor stars in the spheroid 
(Pop II) of our Galaxy, which have metallicities going down to at least 10 -4 and perhaps 
10 -5 of the Solar value [36]. Observations have long shown [37-41] that Li does 
not vary significantly in Pop II stars with metallicities ^ 1/30 of Solar — the "Spite 
plateau" [37]. Recent precision data suggest a small but significant correlation between 
Li and Fe [38] which can be understood as the result of Li production from Galactic 
cosmic rays [39]. Extrapolating to zero metallicity one arrives at a primordial value 
Li/H| p = (1.23 ± 0.06) x 10 -10 [40]. One systematic error stems from the differences 
in techniques to determine the physical parameters (e.g., the temperature) of the stellar 
atmosphere in which the Li absorption line is formed. Alternative analyses, using methods 
that give systematically higher temperatures, yield Li/H| p = (2.19 ± 0.28) x 10 -10 [41] 
and Li/H| p = (2.34 ± 0.32) x 10 -10 [42]; the difference with [40] indicates a systematic 



D/H = (2.78 ±0.29) x 10 



-5 



(1.2) 



Y p = 0.249 ± 0.009. 



(1.3) 
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uncertainty of a factor of ~ 2. Moreover it is possible that the Li in Pop II stars has been 
partially destroyed, due to mixing of the outer layers with the hotter interior [43] . Such 
processes can be constrained by the absence of significant scatter in Li-Fe [38], and by 
observations of the fragile isotope 6 Li [39]. Nevertheless, depletion by a factor as large 
as ~ 1.6 remains allowed [40,43]. Including these systematics, we estimate a primordial 
Li range, as shown in Fig. 1.1: 



Finally, we turn to He. Here, the only observations available are in the Solar 
system and (high-metallicity) H II regions in our Galaxy [44]. This makes inference 
of the primordial abundance difficult, a problem compounded by the fact that stellar 
nucleosynthesis models for 3 He are in conflict with observations [45]. Consequently, it 
is no longer appropriate to use 3 He as a cosmological probe; instead, one might hope to 
turn the problem around and constrain stellar astrophysics using the predicted primordial 
3 He abundance [46]. 

1.3. Concordance, Dark Matter, and the CMB 

We now use the observed light element abundances to assess the theory. We first 
consider standard BBN, which is based on Standard Model physics alone, so N v = 3 and 
the only free parameter is the baryon-to-photon ratio rj. (The implications of BBN for 
physics beyond the Standard Model will be considered below, §4). Thus, any abundance 
measurement determines ?y, while additional measurements overconstrain the theory and 
thereby provide a consistency check. 

First we note that the overlap in the rj ranges spanned by the larger boxes in Fig. 1.1 
indicates overall concordance. More quantitatively, when we account for theoretical 
uncertainties as well as the statistical and systematic errors in observations, there is 
acceptable agreement among the abundances when 



However the agreement is much less satisfactory if we use only the quoted statistical 
errors in the observations. In particular, as seen in Fig. 1.1, D and 4 He are consistent 
with each other but favour a value of r\ which is higher by ~ 2a from that indicated by 
the 7 Li abundance. Additional studies are required to clarify if this discrepancy is real. 

Even so, the overall concordance is remarkable: using well-established microphysics we 
have extrapolated back to an age of ~ 1 s to correctly predict light element abundances 
spanning 9 orders of magnitude. This is a major success for the standard cosmology, and 
inspires confidence in extrapolation back to still earlier times. 

This concordance provides a measure of the baryon content 




(1.4) 



4.7 < 77 < 6.5 (95% CL) . 



(1.5) 



0.017 < n B h 2 < 0.024 (95% CL) , 



(1.6) 
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a result that plays a key role in our understanding of the matter budget of the universe. 
First we note that Ob <C 1, i.e., baryons cannot close the universe [47]. Furthermore, 
the cosmic density of (optically) luminous matter is Oj um ~ 0.0024/i _1 [48], so that 
fig ^> Oi um : most baryons are optically dark, probably in the form of a ~ 10 6 K X-ray 
emitting intergalactic medium [49]. Finally, given that fijyi ~ 0.3 (see Dark Matter and 
Cosmological Parameters reviews), we infer that most matter in the universe is not only 
dark but also takes some non-baryonic (more precisely, non-nucleonic) form. 

The BBN prediction for the cosmic baryon density can be tested through precision 
observations of CMB temperature fluctuations (see Cosmic Microwave Background 
review). One can determine rj from the amplitudes of the acoustic peaks in the CMB 
angular power spectrum [50], making it possible to compare two measures of rj using 
very different physics, at two widely separated epochs. In the standard cosmology, there 
is no change in rj between BBN and CMB decoupling, thus, a comparison of ?7bbn an d 
VCMB is a key test. Agreement would endorse the standard picture while disagreement 
could point to new physics during/between the BBN and CMB epochs. 

The release of the first-year WMAP results was a landmark event in this test of BBN. 
As with other cosmological parameter determinations from CMB data, the derived ?7cmb 
depends on the adopted priors [51], in particular the form assumed for the power 
spectrum of primordial density fluctuations. If this is taken to be a scale-free power-law, 
the WMAP data implies fl B h 2 = 0.024 ± 0.001 or rj 10 = 6.58 ± 0.27, while allowing for a 
"running" spectral index lowers the value to fig/i = 0.0224 ± 0.0009 or rjiQ = 6.14 ± 0.25 
[52]; this latter range appears in Fig. 1.1. Other assumptions for the shape of the power 
spectrum can lead to baryon densities as low as fl^h 2 = 0.019 [53]. Thus outstanding 
uncertainties regarding priors are a source of systematic error which presently exceeds the 
statistical error in the prediction for t]. 

It is remarkable that the CMB estimate of the baryon density is consistent with the 
BBN range quoted in Eq. (1.6), and is in fact in very good agreement with the value 
inferred from recent high-redshift D/H measurements [26] and 4 He determinations; 
together these observations span diverse envirnoments from redshifts z = 1000 to the 
present. However 7 Li is at best marginally consistent with the CMB (and with D) 
given the error budgets we have quoted. The question then becomes more pressing as 
to whether this mismatch come from systematic errors in the observed abundances, 
and/or uncertainties in stellar astrophysics, or whether there might be new physics at 
work. Inhomogeneous nucleosynthesis can alter abundances for a given ??bbn but will 
overproduce 7 Li [54] . Note that entropy generation by some non-standard process could 
have decreased rj between the BBN era and CMB decoupling, however the lack of spectral 
distortions in the CMB rules out any significant energy injection upto a redshift z ~ 10 7 
[55] . Interestingly, the CMB itself offers the promise of measuring 4 He [56] and possibly 
7 Li [57] directly at z ~ 300 - 1000. 

Bearing in mind the importance of priors, the promise of precision determinations 
of the baryon density using the CMB motivates using this value as an input to 
BBN calculations. Within the context of the Standard Model, BBN then becomes a 
zero-parameter theory, and the light element abundances are completely determined to 
within the uncertainties in ??cmb an d the BBN theoretical errors. Comparison with the 
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observed abundances then can be used to test the astrophysics of post-BBN light element 
evolution [58]. Alternatively, one can consider possible physics beyond the Standard 
Model (e.g., which might change the expansion rate during BBN) and then use all of the 
abundances to test such models; this is the subject of our final section. 

1.4. Beyond the Standard Model 

Given the simple physics underlying BBN, it is remarkable that it still provides the 
most effective test for the cosmological viability of ideas concerning physics beyond the 
Standard Model. Although baryogenesis and inflation must have occurred at higher 
temperatures in the early universe, we do not as yet have 'standard models' for these so 
BBN still marks the boundary between the established and the speculative in big bang 
cosmology. It might appear possible to push the boundary back to the quark-hadron 
transition at T ~ ^4qcd or electroweak symmetry breaking at T ~ 1/y/Gp; however 
so far no observable relics of these epochs have been identified, either theoretically or 
observationally. Thus although the Standard Model provides a precise description of 
physics up to the Fermi scale, cosmology cannot be traced in detail before the BBN era. 

Limits on particle physics beyond the Standard Model come mainly from the 
observational bounds on the 4 He abundance. This is proportional to the n/p ratio which 
is determined when the weak-interaction rates fall behind the Hubble expansion rate at 
Tf r ~ 1 MeV. The presence of additional neutrino flavors (or of any other relativistic 
species) at this time increases hence the expansion rate, leading to a larger value of 
Tf r , n/p, and therefore Y p [6,59]. In the Standard Model, the number of relativistic 
particle species at 1 MeV is g* = 5.5 + \N V , where 5.5 accounts for photons and e^, and 
N u is the number of (nearly massless) neutrino flavors (see Big Bang Cosmology review). 
The helium curves in Fig. 1.1 were computed taking N v = 3; the computed abundance 
scales as Z\>bbn — 0.013 AN,, [7]. Clearly the central value for N v from BBN will 
depend on rj, which is independently determined (with weaker sensitivity to N v ) by the 
adopted D or 7 Li abundance. For example, if the best value for the observed primordial 
4 He abundance is 0.249, then, for n\Q ~ 6, the central value for N v is very close to 3. This 
limit depends sensitively on the adopted light element abudances, particularly l^BN- A 
maximum likelihood analysis on n and N u based on the above 4 He and D abundances 
finds the (correlated) 95% CL ranges to be 4.9 < 7710 < 7.1 and 1.8 < N v < 4.5 [60]. On 
the other hand, using the 'low' 4 He and 7 Li gives [61] 7710 < 4.3, and 1.4 < N u < 4.9. 
Similar results were obtained in another study [62] which presented a simpler method 
[10] to extract such bounds based on x statistics, given a set of input abundances. 
Using the CMB determination of n improves the constraints: with a 'low' 4 He, N v = 3 
is barely allowed at 2a [63] but using the 4 He (and D) abundance quoted above gives 
5.66 < 7/10 < 6.58 (n B h 2 = 0.0226 ± 0.0017) and N v = 3.24 ± 1.2 at 95% CL [60]. 

Just as one can use the measured helium abundance to place limits on g* [59], any 
changes in the strong, weak, electromagnetic, or gravitational coupling constants, arising 
e.g. from the dynamics of new dimensions, can be similarly constrained [64]. 

The limits on N v can be translated into limits on other types of particles or particle 
masses that would affect the expansion rate of the Universe during nucleosynthesis. 
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For example consider 'sterile' neutrinos with only right-handed interactions of strength 
Gr < Gp. Such particles would decouple at higher temperature than (left-handed) 
neutrinos, so their number density (oc T 3 ) relative to neutrinos would be reduced 
by any subsequent entropy release, e.g. due to annihilations of massive particles that 
become non-relativistic in between the two decoupling temperatures. Thus (relativistic) 
particles with less than full strength weak interactions contribute less to the energy 
density than particles that remain in equilibrium up to the time of nucleosynthesis 
[65]. If we impose N v < 4 as an illustrative constraint, then the three right-handed 
neutrinos must have a temperature 3(Tj, R /T^ L ) 4 < 1. Since the temperature of the 
decoupled z/r's is determined by entropy conservation (see Big Bang Cosmology review), 
T UR /T Uh = [(43/4)/ i g^(T c j)] 1 / 3 < 0.76, where is the decoupling temperature of the 
i^r's. This requires g*(T^) > 24 so decoupling must have occurred at > 140 MeV. 
The decoupling temperature is related to Gr through (Gr/Gf) 2 ~ (^d/3MeV) -3 , where 
3 MeV is the decoupling temperature for v^s. This yields a limit Gr ^ 10 _2 Gp. The 
above argument sets lower limits on the masses of new Z' gauge bosons in superstring 
models [66] or in extended technicolour models [67] to which such right-handed neutrinos 
would be coupled. Similarly a Dirac magnetic moment for neutrinos, which would allow 
the right-handed states to be produced through scattering and thus increase can be 
significantly constrained [68], as can any new interactions for neutrinos which have 
a similar effect [69]. Right-handed states can be populated directly by helicity-flip 
scattering if the neutrino mass is large enough and this can be used to used to infer e.g. 
a bound of m VT ^ 1 MeV taking N v < 4 [70]. If there is mixing between active and 
sterile neutrinos then the effect on BBN is more complicated [71] . 

The limit on the expansion rate during BBN can also be translated into bounds on 
the mass/lifetime of particles which are non-relativistic during BBN resulting in an even 
faster speed-up rate; the subsequent decays of such particles will typically also change 
the entropy, leading to further constraints [72]. Even more stringent constraints come 
from requiring that the synthesized light element abundances are not excessively altered 
through photodissociations by the electromagnetic cascades triggered by the decays 
[73-75], or by the effects of hadrons in the cascades [76]. Such arguments have been 
applied to e.g. rule out a MeV mass v T which decays during nucleosynthesis [77]; even 
if the decays are to non-interacting particles (and light neutrinos), bounds can be derived 
from considering their effects on BBN [78]. 

Such arguments have proved very effective in constraining supersymmetry. For example 
if the gravitino is light and contributes to the illustrative BBN limit N v < 4 requires 
its mass to exceed ~ 1 eV [79]. In models where supersymmetry breaking is gravity 
mediated, the gravitino mass is usually much higher, of order the electroweak scale; such 
gravitinos would be unstable and decay after BBN. The constraints on unstable particles 
discussed above imply stringent bounds on the allowed abundance of such particles, which 
in turn impose powerful constraints on supersymmetric inflationary cosmology [74-76]. 
These can be evaded only if the gravitino is massive enough to decay before BBN, i.e. 
m 3/2 ~ ^O TeV [80] which would be unnatural, or if it is in fact the LSP and thus stable 
[74,81]. Similar constraints apply to moduli - - very weakly coupled fields in string 
theory which obtain an electroweak- scale mass from supersymmetry breaking [82]. 
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Finally we mention that BBN places powerful constraints on the recently suggested 
possibility that there are new large dimensions in nature, perhaps enabling the scale of 
quantum gravity to be as low as the electroweak scale [83]. Thus Standard Model fields 
may be localized on a 'brane' while gravity alone propagates in the 'bulk'. It has been 
further noted that the new dimensions may be non-compact, even infinite [84] and the 
cosmology of such models has attracted considerable attention. The expansion rate in the 
early universe can be significantly modified so BBN is able to set interesting constraints 
on such possibilities [85]. 
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